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ABSTRACT 

High resolution spectral studies were undertaken at orbital phases (0) 0, 0.25 
and 0.5 on the high-mass X-ray binary (HMXB) Vela X-l using archival Chandra 
data. We present (a) the first detailed analysis of the multiple strong narrow 
emission lines present in 0.5 (b) an analysis of the absorption of the continuum 
in <p 0.5, and (c) the first detection of narrow emission and absorption lines in 
0.25. Multiple fluorescent and H-and He-like emission lines in the band 1.6 - 20 
Angstrom (A) in eclipse are partially obscured at 0.25 by the X-ray continuum. 
The 0.25 spectrum displays 3 triplets, 2 with a blue-shifted resonance (r) line 
in absorption and the intercombination (i) and forbidden (f) lines in emission, 
and shows in absorption other blue-shifted lines seen in emission in eclipse. At 
0.5 the soft X-ray continuum diminishes revealing an "eclipse-like" spectrum, 
however line flux values are around 13-fold those in eclipse. We conclude the 
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narrow emission lines in Vela X-l become apparent when the continuum is blocked 
from line of sight, either by eclipse or by scattering and/or absorption from a wake 
or cloud. The H-and He-like lines arise in warm photoionised regions in the stellar 
wind, while the fluorescent lines (including a Ni Ka line) are produced in cooler 
clumps of gas outside these regions. Absorption of the 5-13 A continuum at 0.5 
may be caused by an accretion wake comprised of dense stagnant photoionized 
plasma inside a Stromgren zone. Multiple fluorescent emission lines may be a 
common feature of the supergiant category of HMXBs. 

Subject headings: X-rays:binaries — pulsars(individual, Vela X-l) — supergiants 
— circumstellar matter — line formation — stars:winds 



1. Introduction 

High-mass X-ray binaries (HMXBs) including Vela X-l have served to illuminate the 
evolution of giant binary stars, to probe their radiation-driven winds, and to study accretion 
processes and the fundamental parameters of compact objects. The 130 known HMXBs fall 
into 2 categories, Be/X-ray and supergiant systems (Liu et al. 2000). Vela X-l is one of 
11 supergiant systems (which include GX 301-2 and 4U 1700-37), while the vast majority 
(around 80%) of HMXBs are Be/X-ray systems detected as transients (Kaper 2001; White, 
Nagase & Palmar 1995). Supergiant systems have a relatively short life of 10,000 years before 
the secondary star explodes in a supernova leaving a neutron star binary or 2 single neutron 
stars (Kaper 2001). Once a massive main-sequence star evolves to a supergiant the Roche 
lobe may become filled and accretion flow becomes sufficient to power a strong X-ray source. 

Vela X-l is an eclipsing system at an estimated distance of 1.9 kpc (Sadakane et al. 1985) 
containing an B0.5Ib supergiant (mass 23 M and radius 34 R & ) and a wind-fed pulsar with 
a period of 283 s (Nagase 1989). The neutron star (NS) orbits in a period of 8.96 days very 
close to the non-synchronously rotating surface of the supergiant with a semi-major axis 
53.6 i? (Barziv et al. 2001), so it is fully embedded in the acceleration zone of the stellar 
wind. The NS mass of 1.9 M Q is the highest known (Barziv et al. 2001; van Kerkwijk et al. 
1995) and the system has been studied at infrared (Huthoff & Kaper 2002), UV (Sadakane 
et al. 1985; van Loon, Kaper & Hammerschlag-Hensberge 2001), gamma (Raubenheimer et 
al. 1994) and at optical (van Kerkwijk et al. 1995; Kaper et al. 1994) wavelengths. It has 
not been detected at radio frequency (detection limit at mean cm radio flux density (mJy) 
< 0.2 (Fender & Hendry 2000). 

Orbital modulation of the spectra have been noted in the UV, X-ray, and optical wave- 
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lengths. Sadakane et al. (1985) noted absorption of Al III and Fe III in the UV, limited to the 
orbital phase (0) 0.5 and higher, and postulated the presence of a trailing wake. Bessell, Vi- 
dal, & Wikramsinghe (1975) reported an absorption component in the Ha spectrum around 
0.6. Kaper et al. (1994) postulated that an absorption component in the optical spectrum 
resulted from a photoionization wake that trails the X-ray source, rather than the accretion 
wake which does not sufficiently obscure the supergiant. Smith (2001) studied an asymmet- 
ric dip in the UV lightcurve at 0.46-0.7 and concluded spectral ratios (predip/dip) contain 
information about the mechanism responsible for the variability. Absorption of soft X-rays 
during the 0.5 and higher was first reported in Tenma data (Nagase 1989). Feldmeier et al. 
(1996) showed that a trailing photoionization wake can explain the asymmetry in the X-ray 
band during the eclipse ingress and egress. 

Sako et al. (1999) using ASCA data demonstrated that the spectrum in eclipse was 
comprised of recombination lines and radiative combination continua (RRC) produced by 
photoionization in the stellar wind; they suggested the fluorescent K-shell lines from near- 
neutral atoms indicated that the X-ray -irradiated portion of the wind in Vela X-l consists 
of cool dense clumps in a hotter, more ionized gas. 

Schulz et al. (2002a) presented the first high-resolution spectrum during eclipse, and 
using RRC (including Ne X RRC) deduced a temperature of 1.2 x 10 5 K for the plasma. The 
spectrum showed clear evidence of photoionization processes, and the resonance (r) lines in 
the He-like triplets were of roughly equal strength to the forbidden (f) lines, consistent with 
resonant scattering. 

We unexpectedly found that in 0.5 the emission lines were stronger than in eclipse by 
an order of magnitude, associated with the near disappearance of the soft X-ray continuum. 
This result, together with the unusual finding of multiple fluorescent emission lines, led us 
to further investigate emission line production in HMXBs by a study of Vela X-l comparing 
eclipse, 0.25 and 0.5 using high-resolution Chandra datasets. We outline the data 
reduction proceedure in section 2, give our results in section 3 and present the discussion in 
section 4. 



2. Observations 

An analysis of Vela X-l in eclipse (Chandra Dataset Obsld 102) has previously been 
reported (Schulz et al. 2002a). Three newer datasets obtained during one binary orbit (Figure 
1) in 2001 with the Chandra High Energy Transmission Grating Spectrometer (HETGS) 
representing the orbital phases eclipse (0 0), 0.25 and 0.5 were available from the public 
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archive (ObsIDs 1926, 1928, 1927, respectively (see Table 1). These datasets have not been 
previously reported, except for a specialized study of pulse-phased resolved emission of the 
Fe Ka line (Paul et al. 2002), and a summary report in a review of stellar winds in 5 HMXBs 
(Sako et al. 2003). The telescope operated in "faint" data mode, with a timed readout (frame 
time) of 1.74 seconds. The HETGS employs two sets of transmission gratings: the Medium 
Energy Grating (MEG) with a range of 2.5-31 Angstroms (A) (0.4-5.0 keV) and the High 
Energy Grating (HEG) with a range of 1.2-15 A (0.8-10.0 keV). The resolution of the first- 
order HEG spectrum is AA = 0.012 A, and MEG spectrum AA = 0.023 A. The dispersed 
spectra were detected with the ACIS-S linear array of 6 CCD detectors. The HETGS/ACIS- 
S combination provides an undispersed (zeroth order) image and dispersed spectra from the 
gratings. The zeroth-order image of the non-eclipse datasets were affected by photon pileup 
and were not used in spectral analysis. 

A complete re-processing was undertaken using Ciao software (version 2.3 with Caldb 
version 2.22) (including the tools acis_process_events, tgdetect etc) (Chandra X-ray Cen- 
ter Data analysis 1 ) to produce grating ancillary response files (ARFs) and redistribution 
matrix files (RMFs). For spectral response the files acisheglD1999-07-22rmfN0004.fits and 
acismeglD1999-07-22rmfN0004.fits were used for the HEG and MEG respectively. 

Measurement of emission lines in photoionized spectra involves consideration of the 
contribution to the lines by a continuum, and using Isis (version 1.0.50) lines were fitted 
with a gaussian plus a polynomial to fit the narrow regions adjacent to the line. 

Initial inspection of the 0.25 spectrum showed a jagged appearance with few apparent 
narrow features apart from the Fe Ka line. For all datasets a fast fourier transform routine in 
Isis was used to smooth the spectral curve, which was viewed "in counts" to assess statistical 
quality and find instrumental features, "in flux density" to compare HEG and MEG spectra, 
and in "delta-chi" to see the relative significance of deviations. This procedure identified a 
number of previously unreported emission and absorption features in the 0.25 spectrum, 
including 3 He-like triplets. " Features" that do not appear both in HEG and MEG spectra 
are likely noise, or are more prominent in the HEG spectrum due to its higher resolution 
(Figures 2,3). Broad systematic flux differences are due to residual calibration errors. For 
0.25 spectra, the statistical uncertainties are about 6% and 11% for MEG and HEG, 
respectively, near 7.5 A, and rise to about 18% and 32% near 11.5 A. 

Pileup in a Charge-Coupled Device (CCD) occurs when there is coincidence of two or 
more photons per CCD frame-time within an event detection cell. Two or more photons 
interact in a detection cell of the detector resulting in a pulse height that is roughly the 
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sum of the pulse heights of the individual photon events. This causes a lower event detection 
rate and distortion of the observed spectrum towards higher energies (Davis 2001). The 
zero-order images of ObsID 1928 and 1927 showed visual evidence of pile-up. We checked 
for pileup in Vela X-l (ObsIdl928) first dispersion order very carefully, as its flux and mean 
MEG count rates are comparable to another X-ray binary 4U 1543-624 (Obsld 702) which is 
known to have mild pileup in the first order. The important metric in determining pileup is 
the counts/frame/cell, where a cell is a 3 pixel bin along the spectral dispersion axis (Dewey 
D 2 ) 

Vela X-l has 0.0345 and 4U 1543-624 has 0.045 counts/frame/cell, that is 4U 1543-624 
has counts/frame/cell 30% higher than Vela X-l. Dewey argues the value of this metric indi- 
cates the magnitude of pileup effects: for values of 0.01 and less pileup effects are negligible; 
for values of 0.01 to 0.1 pileup may produce changes in event rates less than 10%; for values 
above 0.1 pileup can be severe and the data quality dubious. 

The lack of pileup in Vela X-l is demonstrated by a careful comparison of the HEG and 
MEG datasets in the range 1.7-15.2 A. There is close agreement in the first-order spectra. In 
a graded series of 5 spectra with varying degrees of pileup, all spectra with pileup (including 
4U 1543-624, the mildest case of the 5) show lack of close agreement between the MEG and 
HEG (J Davis 3 ). The MEG has a larger effective area than HEG, and in the event of pileup 
the MEG will be piled (or more piled) than HEG, and the MEG spectrum saturates and 
does not follow (shows lower flux) than the HEG. 

The assessment of absorption features in the 0.25 spectrum found photoelectric neutral 
K absorption edges from S, Si and Mg in cold gas. An optical depth (tau) of the neutral 
edge was assessed: 

T — ^(fhigh/ flow)', 

where " fhigh" is the X-ray flux on top of the edge at the high wavelength side and 
"flow" is the X-ray flux at the bottom of the edge at the low wavelength side (Schulz et al. 
2002b). This optical depth is directly proportional to neutral absorption column densities 
and assesses the total amount of photo-electric absorption, which includes absorption in the 
vicinity of the source as well as the interstellar medium in the line of sight. 

Smith (2001) in studying the dip in the UV lightcurve at 0.46-0.7 compared spectra 
predip to dip to obtain a "ratioed spectrum": 



2 http://space. mit.edu/HETG/technotes/pileup/pileup_9912. html 

3 http://cxc. harvard.edu/ccw/proceedings/02_proc/presentations/j_davis/pileup. html 
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Fobs (A) — 1 — fmin(^)/fmax(^) 

Between 5 and 13 A, the strong X-ray continuum noted in 0.25 was markedly at- 
tenuated in 0.5. We consider a simple hypothesis in which an absorber (accretion wake 
associated with Bondi-Hoyle accretion) causes absorption at 0.5, but is not in the line-of 
sight at 0.25 (Figure 4). We estimated the "ratioed spectrum," and also modeled an optical 
depth of the absorption using: 

T = ln(f phase0 .25/ fphaseO. 5)- 

3. Results 

The 3 datasets were observed during 1 orbit of the companion by the NS (Figure 1) 
and no flaring behaviour occurred. The lightcurves from the 3 observations are compared in 
Figure 5, in which pulsed emission (period 283 s) may be discerned. 

The 0.5-10 keV spectra were fit with a model of 2 absorbed powerlaws, a scattered and 
direct component (Sako et al. 1999). The photon index was fixed to 1.7 in both powerlaws 
in eclipse, with the second powerlaw index allowed to vary in the other phases (see Table 1, 
Figure 6). The total flux in 0.5 is reduced compared to 0.25 (Table 1, Figure 2). 

Details of selected spectral lines and RRC are set out in Tables 2, 3,& 4. A rich mosaic 
of 10 H-like, 12 He-like and 15 fluorescent lines (including the unusual Ni Ka line), and 2 
RRC were detected in the eclipse phase. Our results in eclipse phase are in general agreement 
with those reported by Schulz et al. (2002a), and were not significantly different from ObsID 
102. 

Comparisons of 2 spectra (from 0.25 and 0.5) are presented in Figure 2 (wavelength 
range 1.5-14.3 A) and Figure 3 (1.6-2.0 A). 

The 0.25 displays a mix of narrow features in emission and in absorption (Tables 2, 3 
; Figures 2,3), and 5 photo-electric absorption edges were detected (Table 5). The narrow 
emission and absorption lines in this phase (including three He-like triplets) have not been 
previously reported. Two triplets have blue-shifted r lines seen in absorption, while the 
intercombination (i) and f lines are seen in emission (Figures 2, Table 3). 

Comparing 0.5 to eclipse, the H- and He-like emission lines are 13 fold stronger at 
0.5 (mean value for comparisons on 11 lines where firm confidence intervals established). 
The Fe Ka is 21 fold and Fe Kf3 31 fold stronger. The Ca and Ar fluorescent lines seen in 
eclipse are not found in 0.5, probably due to the stronger direct continuum below 5 A in 



-7- 



0.5. The Ni Ka line and a Fe K edge are shown in Figure 3. Vela X-l is the only XRB 
with a Ni Ka line reported, however our study of GX 301-2 (in progress) has detected this 
line (Table 2). 

In and 0.5 the r and f lines of the He-like triplets have roughly equal strength. 
However in 0.25 the triplet r line is a little weaker than the f line (Ne IX) or is in absorption 
(MgXI, Si XIII) (Figure 2). 

The He and H-like emission lines showed a red-shift at (of around 300 km s -1 ), 
and at 0.25 (of around 200 km s _1 ). However lines seen in absorption at 0.25 have a 
blue shift of around 300 km s" 1 , consistent with the orbital velocity of the NS. At 0.5 the 
emission lines have a blue-shift (of around 150 km s" 1 ). 

The dramatic attenuation of the X-ray continuum at 0.5 compared to 0.25 is modeled 
both as a spectral ratio (0 0.25/0 0.5) and as an optical depth in Figure 7. The optical depth 
increases steeply between 3 and 6 A, then gradually declines from 6-13 A. 

4. Discussion 

These results provide insights into how narrow lines are produced in HMXBs, and 
particularly into multiple fluorescent emission and orbital changes in emission line strength. 
The main features to be discussed below are: 

(a) changes with orbital phase in emission lines and the strong narrow emission lines 
present in 0.5, the first detection of narrow emission and absorption lines in 0.25, and 
an analysis and measurement of the absorption of the X-ray continuum in 0.5. 

(b) new aspects of the multiple fluorescent emission in Vela X-l, including the identi- 
fications of the Ni Ka line in Vela X-l, and also in GX301-2; the variation in line strength 
of Fe and Ni Ka with orbital phase, attributed (following Matt et al 1997) to an accreting 
system where the direct radiation is blocked, and only reflected radiation is visible; and the 
finding that multiple fluorescent emission lines may be a common feature of the supergiant 
category of HMXBs. 

(c) Doppler shifts of emission lines in 0.5 indicate a slowing of the wind associated 
with the presence of the X-ray source and formation of the Stromgren zone and accretion 
wake. Two triplets in 0.25 show blue-shifted resonance lines seen in absorption while 
the intercombination and forbidden lines are seen in emission, with the blueshifted line in 
absorption having a different source to the lines seen in emission. 
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The X-ray luminosity of Vela X-l is known to demonstrate considerable variation on 
time intervals from hours to days (Kreykenbohm et al. 1999; van Kerkwijk et al. 1995); there 
may be pulse-to-pulse variations even though the pulse profile averaged over many pulses is 
quite stable. Clearly our Chandra results (where flux change, and absorption change between 
phases are carefully measured) relate to only the observed orbit. 

4.1. Photoionized Plasma 

In Vela X-l spectral features indicating plasmas in photoionization equilibrium (PIE) 
include (a) multiple fluorescent ions (b) narrow RRCs (Ne X and Ne IX) (c) relatively weak 
or apparent absence of iron L-shell emission in presence of K-shell emission from lighter 
elements (S, Si, Mg, Ne and O) (Liedahl et al. 2001), and (d) He-like triplet ratios that are 
consistent with photoionization and resonant scattering (involving photoexcitation) in a PI 
plasma. For Ne X the Lyman series (Lya:P:^:8) have line ratios of 1.0:0.35:0.17:0.12 in eclipse 
and 1.0:0.46:0.34:0.13 in 0.5. These flux ratios support the existence of photoexcitation 
since the fluxes of the higher lines in the series are higher than expected from a purely 
recombining or purely collisional plasma (Schulz et al. 2002a). 

PIE are characterised by processes of photoionization (absorbed photon ejects electron 
from ion) and recombination (electron recombines with ion and emits photon or excites an- 
other electron). In PIE the ratio of the X-ray flux to the electron density is sufficiently 
high that collisional ionization processes are negligible relative to photoionization processes 
(Heeter et al. 1999). Spectra of PIE consist of either recombination emission, or fluorescent 
emission, or a composite of these (Liedahl 1999). The dominance of fluorescence over re- 
combination, or vice versa, is related to the charge state distribution, with highly ionized 
plasmas being recombination-dominated, while colder, less ionized plasma are fluorescence- 
dominated. 

The He-like triplets have r and f line fluxes of roughly equal strength (except in the <j) 
0.25, see section 4.4) (Table 3, Figure 2). This is consistent with resonant scattering, defined 
as photoexcitation and subsequent decay by Wojdowski et al. (2003). The strong r line 
in eclipse in a photoionized plasma occurs because the Lja lines of H-like ions and the 
n = 2 — >■ 1 r lines of He-like triplets have large oscillator strengths (unlike the i and the f 
lines). Line emission of lines with large oscillator strengths is enhanced in the wind of an 
X-ray binary by photoexcitation from radiation from the compact source. However outside 
of eclipse photons from the NS may be scattered out of the line of sight to the observer, 
reducing the line flux of the r line of the He-like triplet (and the Lya lines), and if sufficient 
photons are scattered the r line may be seen in absorption, as noted in 0.25 with Si XIII and 
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Mg XI (Table 3, Figure 2). In 0.5 the absorption of the continuum creates an "eclipse-like" 
spectrum and triplets in and 0.5 display similar ratios. 

Narrow RRCs from Ne X and Ne IX (at 9.10 A and 10.37 A) were detected in eclipse and 
0.5 and indicate a plasma temperature of 1.2 xl0 5 fT (Table 4). The measured Lya/RRC 
flux ratio for Ne X was 1.9 in eclipse and 2.1 in 0.5, which is not consistent with the ratio 
of 1.3 predicted for a purely recombining plasma, and implies photoexcitation is present 
(Schulz et al. 2002a). 

4.2. The fluorescent lines. 

There are 9 NS systems in the OB-supergiant category of HMXBs, to which Vela X-l 
belongs, and 2 black hole systems (Kaper 2001). Of 5 systems for which high- resolution X-ray 
spectra are available, multiple fluorescent lines are present in 3: Vela X-l (this study, Schulz 
et al. (2002a)), GX 301-2 (our study in progress) and 4U 1700-37 (Boroson et al. 2003). Our 
study of GX 301-2 reveals Ca II- VII at 3.358 A, Ar VI- IX at 4.18 A, Si IX at 5.316 A, 
S V at 3.367 A and Si V at 7.116 A, in addition to Ni and Fe Ka lines. Fe Ka transitions 
are common in photoionized X-ray spectra, but otherwise a multitude of fluorescent lines 
is unusual in X-ray astronomy. Matt, Fabian & Reynolds (1997) suggest that fluorescent 
emission from elements other than iron from neutral, externally X-ray illuminated matter 
may occur if the direct radiation is blocked and only reflected radiation is visible; then 
neutral ions of Ne, Mg, Si, S, Ar, Ca, Cr and Ni may be detected, and the line flux increases 
with increasing inclination angle of the gas slab. Thus in both Vela X-l and GX301-2 the 
Fe Ka flux may vary in orbit as the angle changes, and when it is maximal the Ni K a line 
(and other fluorescent lines) are most apparent (Table 2). 

Photoionization and Compton ionization by high-energy X-ray or Gamma-ray photons 
are effective ways to remove inner shell electrons from an atom or ion, and the vacancy in 
the inner shell is filled by electron cascades from higher level shells of the ion. The process 
involves either a radiative transition - fluorescence - or a non-radiative transfer of energy to 
an electron (Auger effect), so the fluorescence yield, defined as the probability that radiation 
is emitted during the filling of the vacancy, is less than 1. The probability scales as Z 4 , so 
high-Z Ka lines such as Fe are more common (Liedahl et al. 2001). The Ni Ka emission 
line in Vela X-l was first detected by Sako et al. (1999), and is present in the scattered 
component of the eclipse and 0.5 spectra. This line is rarely detected in X-ray binaries 
or X-ray astronomy in general, and the first unambiguous detection in an active galactic 
nucleus was reported in 2003 (Molendi, Bianchi & Matt 2003). 
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A Ka line occurs where an inner-shell ionization event causes a vacancy in the Is 
shell, and the electron making a transition to fill the vacancy is a 2p electron. By contrast 
in collisional plasmas the processes of collisional excitation and excitation-auto-ionization 
mostly affect the outer electrons, and are less often associated with inner-shell auto-ionization 
(Kaastra & Mewe 1993). 

Fluorescent lines in Vela X-l are an indicator of cooler clumps in the wind. Sako et 
al. (1999) developed a differential emission measure distribution to model Vela X-l, which 
has well-constrained system parameters, and achieved a realistic reproduction of the ASCA 
spectrum. The model had to accomodate lines from near-neutral Mg, Si, S, Ar, Ca, Fe, 
and Ni. In a simple homogenous wind model the local ionization parameter in the emission 
region was far too high to support near-neutral ions. The fluorescent lines must be produced 
in dense clumps or filaments in the wind that may be 100 times denser than in the ambient 
wind, and have local ionization parameters that are up to 2 orders of magnitude lower so they 
can support ions of low charge state. Liedahl (1999) has described in irradiated accretion 
disks how hard X-ray photons can "burrow" into cool material, and are absorbed through 
K-shell photoionization producing fluorescent emission. 



4.3. Doppler Shifts of Emission lines 

The orbital velocity of the NS was calculated at 300 km s _1 , consistent with Sako et al. 
(1999). In 0.25, the lines seen in absorption (directly in line of sight of the NS) generally 
had a blueshift consistent with this value (Tables 2, 3). However other lines seen in emission 
(as a result of scattering) had a redshift of the order of 200 km s -1 , indicating their source 
in ions in part of the stellar wind with a relative bulk motion of around 500 km s _1 . 

In eclipse the H- and He-like ions demonstrate a modest redshift, and in <fi 0.5 they are 
blueshifted. The lines are produced in a process (resonant scattering) whereby ions in the 
wind absorb a photon from the X-ray source, and emit a photon in the line of sight to the 
observer. In eclipse these ions are located in regions of the wind (from the supergiant) which 
are receeding, and out-of eclipse are approaching the observer. The mean blueshift was 179 
km s _1 in 13 lines in <fi 0.5 where fits with narrow (< ±200 km s _1 ) confidence limits were 
obtained. This indicates that the supergiant's wind velocity is less than expected, since at 
the distance of the NS from the supergiant's surface, it should be around 860 km s" 1 , that 
is around 50% of the wind terminal velocity of 1700 km _1 (Dupree et al. 1980)). The result 
suggests that the high photoionization conditions effectively quench the driving force and 
acceleration of the wind, contributing to the formation of a Stromgren zone (see section 4.5 
below). 
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Precise line-shifts were difficult to obtain in the fluorescent ions, with the exception of 
the Fe Ka line. The Fe Ka line does not follow the same pattern as H- and He-like ions: 
although there is a small redshift apparent in eclipse, in 0.5 there is no blueshift as found 
in the H- and He-like ions. This is consistent with the above argument that the fluorescent 
ions arise in different regions of the stellar wind. The Fe Ka line has a redshift around 200 
km s~ x in and 0.5, and a small blueshift at 0.25. This shift is consistent with lightly 
ionized Fe VIII- X (House 1969), and the orbital motion of 300 km s _1 may explain the 
blueshift at 0.25. 



4.4. Absorption and Emission lines at 0.25 

Haberl, White & Kallman (1989) noted that when the X-ray source of HMXBs is viewed 
through the stellar wind captured by the compact object, absorption K-edges such as O, Si, 
S, and K- and L- edges of Fe are seen. Many of these edges are apparent in 0.25, but only 
the K-edges of S and Fe are seen in 0.5 (Table 5). 

Several lines seen in emission in eclipse appear in absorption at 0.25 (Ne X Ly7, the 
r lines of the Mg XI (Figure 2, Table 3) and Si XIII triplets (Table 3), and Mg V (Table 2). 
These features only became clear when the smoothing technique was applied. Smoothing (or 
binning) on the scale of the resolution make weak absorption and emission features easily 
visible. The photons from the X-ray source in the line of sight may be scattered or absorbed 
by the stellar wind en route to the observer. Miller et al. (2002) also noted these triplets 
with r in absorption and i and f in emission in another supergiant HMXB, the black hole 
system Cyg X-l. We expect the phenomenon may also be found in other X-ray binaries such 
as 4U 1700-37, when spectra at various orbital phases are available. Some emission lines are 
seen at 0.25 as a result of resonant scattering (Ne X and Mg XII Lja and Ly/3 lines, Ne IX 
triplet), although others may be obscured by the continuum. The Doppler shifts of the lines 
seen in absorption and those seen in emission are different reflecting their different regions 
of origin. 

In 0.5 no absorption features are detected, apart from edges at 5.0 and 1.74 A(Figure 
3). Phase 0.25 has numerous absorption lines and edges, so a further study between these 2 
phases may be warranted to find a transition phase with more absorption features. 
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4.5. Strong Emission Lines at 0.5 

It is generally considered in HMXBs such as Cen X-3 and Vela X-l that the emission 
lines are most prominent in eclipse, but persist outside of the eclipse phase of the orbit 
(Miller et al. 2002). We consider the marked absorption of the continuum at 0.5 (Figure 
7) to be responsible for the unexpectedly strong emission lines seen in this phase. In eclipse 
only X-rays which have been reprocessed in a visible part of the wind are observed. The 
0.5 line emission displays the same features as in 0, only fluxes are stronger. Thus at 
0.5 there is an "eclipse-like" condition with the disappearance of the X-ray continuum in the 
band 4-12 A associated with the detection of the strong line emission due to scattering of 
photons in the stellar wind. 

SMC X-l is a further example of stronger line emission apparent when the continuum is 
blocked out of eclipse (Vrtilek et al. 2001). In SMC X-l when out-of-eclipse in the high state 
the X-ray source is not blocked and the observer sees only a continuum (similar to the Vela 
X-l 0.25 spectrum), but when the source is blocked by precession of the disk (signified by 
the low-state) the line emission is apparent (Vrtilek et al. 2001). Vela X-l has no precessing 
disk, but the absorption of soft X-rays and other radiation bands around 0.5 has been 
reported by many investigators. 

Liedahl et al. (2001) note that if local material intervenes with our line of sight to the 
continuum source, then absorption may produce observable photoelectric edges or otherwise 
alter the shape of the continuum. The X-ray irradiated gas in the vicinity of the continuum 
source radiates lines and continua characteristic of the local conditions and may intercept 
continuum radiation at energies corresponding to atomic level separations. Thus a photon 
can be absorbed, and subsequently re-emitted at nearly the same energy, but with a new 
direction, so a subtraction from the photon flux propagating toward the observer may occur. 

We postulate that a wind structure in the line of sight accounts for the increased optical 
depth in the 3-12 A band in 0.5. Two distinct structures are described in the various reports 
of absorption (including those that find a higher absorption after 0.5 than before it). These 
are the photoionization wake, and the accretion wake (Figure 4). The X-ray source creates a 
Stromgren zone in which the high photoionization conditions effectively quench the driving 
force and acceleration of the wind. The stagnant flow inside the ionization zone collides 
with the accelerating wind at the trailing border to form a photoionization wake, as shown 
in Figure 4. In the UV investigators have concluded the dip in the lightcurve around 
0.5 is due to the photoionization wake, because the source is large (the supergiant) and the 
accretion wake could not sufficiently occlude it. By contrast the X-ray source is very compact 
and is readily occluded by the accretion wake at 0.5 (Figure 4). Dumm et al. (2000) has 
reported how an accretion wake in an X-ray nova has modified the lightcurve in the UV. 
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Supergiant stars have prodigious mass loss through stellar winds, typically 10~ 6 M Q per 
year. Wind material departs uniformly from the surface of the star, and accretes onto the 
NS if it passes within a critical distance forming an accretion column. Adjacent material 
that is perturbed in its flow but is not accreted forms an accretion wake. 

The discovery of this wind structure in Vela X-l may be helpful in interpreting emission 
lines in other X-ray binaries. 4U 1700-37 is also reported to have strong emission lines 
out-of-eclipse (0 0.7) (Boroson et al. 2003). These authors consider the possibility that the 
emission lines appear because the X-ray source is seen "through an enhanced-density gas 
stream" . However it may be argued that the lines are seen via scattered photons. Analogous 
to the case of Vela X-l, the lines may be seen in emission because the continuum X-rays 
approaching directly from the source are absorbed by an accretion stream or wake, allowing 
the scattered photons to be viewed. Haberl et al. (1989) have modelled the accretion stream 
and wake in 4U 1700-37, which was found to act as an absorber at 0.6. 



4.6. Source of X-rays 

Raubenheimer et al. (1994) has suggested that X-rays and gamma rays from Vela X-l 
are produced by the same mechanism, that is by a beam of charged nucleons accelerated in 
the strong magnetic field at some distance from the surface of the NS. The beam interacts 
with accretion material to form mesons which decay to gamma rays which escape when the 
production region is far enough from the NS so they are not depleted by pair production 
in the strong magnetic field, and the density is low enough to avoid absorption. However, 
Liedahl et al. (2001) and Wojdowski et al. (2003) prefer resonant scattering of X-rays from 
an accretion "hot spot" . This study provides further support for resonant scattering. 

Vela X-l may have an accretion column, as well as an accretion wake. Vela X-l's X-ray 
signal is strongly pulsed, indicating there is a hot spot on the NS surface. A hot spot is 
predicted from the cyclotron scattering resonance features with a cyclotron resonance high- 
energy cutoff of 15-20 keV, corresponding to a magnetic field strength of 2.3 xlO 12 Gauss 
(White et al., 1995). The hot spot is a thermal radiator and the pulse profile is determined by 
the variation in the projected area of the hot spot as the NS rotates (Lyne & Graham-Smith 
pp 131-143, 1998). The pulse has high amplitude (amplitude 45% in the 1.2-2.3 keV band 
and 55% in the 13.5-18.5 keV band of the time-averaged signal), and a notch or "double pulse 
shape" in the main pulse in the X-ray and UV pulse profiles (shown in Figure 7 of Boroson 
et al. (1996)). Pulse profiles at different energies are further determined by shadowing of 
surface radiation by the accretion column and gravitational focusing of column radiation. A 
notch in the main pulse may result where the accretion column is optically thick, causing 
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scattering of radiation propagating vertically through the accretion column. 

Insight as to how these complex structures interact to produce a strong absorption of 
the 5-12 A band may be provided by study of pulse-phase resolved emission. In a study of 
the HEG spectra Paul et al. (2002) demonstrated weak phase dependence of the Fe Ka line 
and suggested that this represented an anisotropic velocity distribution; further studies on 
additional lines are in progress. 

5. Conclusions 

X-ray emission in Vela X-l results from capture and accretion of gas in the stellar wind 
of the supergiant star onto the hot spot on the NS. Wind material departs uniformly from the 
surface of the star, and accretes onto the NS if it passes within a critical distance forming an 
accretion column. Adjacent material that is perturbed in its flow but is not accreted forms 
an accretion wake. The X-rays ionize and heat the surrounding gas, and are reprocessed in 
the stellar wind. Spectra at orbital phases 0, 0.25 and 0.5 reflect the processes of mass loss 
in the stellar wind of the supergiant and mass accretion onto the NS. 

The Doppler shifts at <fi 0.5 indicate slowing of the supergiant's wind velocity in the 
vicinity of the NS, consistent with X-ray photoionization and destruction of ions with strong 
UV resonance transitions that otherwise drive the stellar outflow. The X-ray source creates 
a Stromgren zone in which the high photoionization conditions effectively quench the driving 
force and acceleration of the wind. 

Spectra at <fi 0.25 reveal simultaneous evidence of emission lines produced by scattering, 
and absorption lines that are produced in different regions of the binary system. 

Narrow emission lines in Vela X-l become apparent when the continuum is blocked 
from line of sight, either by eclipse or by scattering and/or absorption from a wake or cloud. 
The X-ray source is very compact and is readily occluded by the accretion wake at <fi 0.5, 
accounting for the absorption of soft continuum. Absorption of the continuum at <p 0.5 
has allowed detection of scattered radiation and emission lines whose strength is an order of 
magnitude higher than those in the eclipse phase. The eclipse and 0.5 spectra have a similiar 
appearance, and the ratios of lines in the He-like triplets are consistent with photoionization 
and resonant scattering (involving photoexcitation) in a PI plasma. 

Multiple fluorescent emission lines have been detected in 3 of 5 supergiant HMXBs for 
which high resolution X-ray spectra are available. This study suggests fluorescent emission 
from externally X-ray illuminated matter may occur from elements other than neutral iron 
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if the direct radiation is blocked and only reflected radiation is visible, when neutral ions of 
Ne, Mg, Si, S, Ar, Ca, Cr and Ni may be detected. 

Differences in Doppler shifts of various lines indicate different origins of lines, from direct 
radiation from the surface of the NS (lines seen in 0.25 in absorption) and from adjacent 
regions in the stellar wind from scattered radiation. 
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Fig. 1. — RXTE ASM lightcurve of 2 orbits of the neutron star in Vela X-l, showing the peri- 
ods (epochs) of the 3 Chandra datasets ( phases (0) 0, 0.25 and 0.5) within one 8.96 day orbit. 

Fig. 2 NEXT PAGE Comparison of smoothed spectra, showing the strong contin- 
uum in phase (0) 0.25 (upper spectrum) that partially obscures the emission lines apparent 
in phase 0.5 (panels 2 and 3). There are 2 curves for each phase, representing the high and 
medium energy grating (HEG and MEG) spectra. In 0.25 HEG (black line) with higher 
resolution has more sharp unresolved features than MEG (red line) , and its lower area gives 
it larger statistical uncertainty. For 0.25 spectra, the statistical uncertainties are about 
6% and 11% for MEG and HEG, respectively, near 7.5 A, and rise to about 18% and 32% 
near 11.5 A. There are three lines (Mg XI triplet r, Ne X Ly7 and Ly/3) in absorption in 
0.25 and in emission in 0.5. The inverted Mg XI triplet r is blueshifted compared to the i 
and f lines. 
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Fig. 3. — (f) 0.5 smoothed spectrum, showing Fe Ka, Fe K/3 and Ni Ka lines, and Fe K edge. 
HEG is green and MEG is red. 



-21 - 



photoionization wake 




phase 0.25 



Fig. 4. — Line of sight at orbital phases 0.25 and 0.5 of the neutron star and its accretion 
wake and photoionization wake. The compact source on the neutron star may be occluded 
by the accretion wake at <ft 0.5. Wind material within a critical radius accretes onto the 
neutron star (curved arrows), while material peturbed in its flow that does not accrete forms 
an accretion wake. 
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Fig. 5. — X-ray lightcurves from the eclipse and (j) 0.25 and </> 0.5 datasets in which pulsed 
emission (period 283 s) is apparent. 
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Fig. 6. — Fitted spectra from 3 phases (eclipse, 0.25 and 0.5), using a best-fit model of 
2 absorbed powerlaws. 
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Fig. 7. — Assessment of continuum absorption in 0.5 by comparison with 0.25, using 
"spectral ratio" and optical depth models (see Observations section for equations). The 
optical depth model (top curve) increases sharply from 2-6 A, then gradually declines 6-13 
A, then declines sharply, while the spectral ratio curve shows a gradual increase from 2-4 A, 
and then is completely flat from 4-18 A. 
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Table 1. Continuum Parameters in Eclipse, 0.25 and 0.5 



Parameter 


Eclipsc(ObsId 102) Eclipsc(ObsID 1926) 


<t> 0.25 


4> 0.5 


Scattered 












0.5 


0.5 


0.5 


0.5 


Powerlaw norm b 


0.00060 


0.00064 


0.0083 


0.0066 


Powcrlaw PI 


1.7 


1.7 


1.7 


1.7 


Direct 










N H 


39.9 


39.0 


1.2 


12.2 


Powerlaw norm 


0.0081 


0.0074 


0.22 


0.12 


Powerlaw PI 


1.7 


1.7 


0.92 


0.60 


Flux(0.5-10 keV) c 


1.17c-ll d 


1.13e-ll 


3.16c-9 


2.04e-9 


Datasct Details 


ObsID 


102 


1926 


1928 


1927 


Start Date 


2000-04-13 


2001-02-11 


2001-02-05 


2001-02-07 


Duration (ks) 

<t> 


30 


85 


35 


30 


0.98-0.02 


0.97-0.08 


0.23-0.27 


0.48-0.52 



a in units of 10 cm 

b in units of 10 -5 photons cm~ 2 s~ 1 

c in units of ergs cm _2 s _1 

d Schulz ct al. (2002a) reported a measurement of 0.88e-ll ergs cm _2 s _1 on this dataset 
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Table 2. Fluorescent emission lines at 3 orbital phases in Vela X-l. Only a few flourescent 

lines are found in 0.25, including the Mg V line which is blueshifted and seen in 
absorption. At the bottom of the table a comparison is made of the 2 fluorescent Fe K lines 
and the rare Ni Ka line in Vela X-l and GX 301-2. 



Ion a Expected A Ob 1926 


Ob 1928 


Ob 1927 


Ratio 


A 

Ca II- VII 3.359 
Ar VI-IX 4.186 
S IX 5.320 
S V 5.367 
Si XI-XII 6.750-6.813 
Si X-XI 6.813-6.882 
Si IX 6.947 
Si VIII 7.007 
Si VII 7.063 
SiV 7.117 
Mg VIII-IX 9.436-9.544 
MgV 9.814 
Mg II 9.890 


<t> 

Shift b Flux c 

-119 0.5 
847 0.2 
198 0.5 
34 0.5 
0.1 
0.3 

Olu -490 u -°-0.36 

-U9 3 _% 8 0.3 
-85 0.2 

-763 0.1 
-217 0.5 


4> 0.25 

Shift Flux d 
63 2 _3 2 ° 51 11.6 

-578i°3 64 -3.lL-!. 8 


(p 0.5 

Shift Flux 

120 2 J| 27 H-5g;| 
-82 3.9 
-377 2.2 
-150 3.8 
1.8 
0.9 

-1028 2 7 1 5 37 3.0 
-396 1.0 
-527! 2 2 1 49 1.6 
2076 14 2 9 33 3.2 

^•ii.O 
539 0.6 

-10 0.6 


Flux 
1927/26 

22.6 

20.5 

4.2 

7.1 

18.0 

3.0 

2.9 

3.2 

8.2 

5.4 
1.3 


Neutral Fe and Ni Fluorescent Lines in Vela X-l 


Ni Ka 1.660 
Fc Kf3 1.756 
Fe Ka 1.937 


-328 4.860 
-47 9 _°| 82 2.2ll 5 
210% 7 18.1 2 _| 6 


41?. 8 3 4 i5 59 2 _ 8 28 
-67%, 221" 19 


1 03 369 74 52 
iu^_ 374 ' 4 -33 

— 126_ 325 52_ 2g 

119i 8 48 374 2 J 27 


15.1 
31.1 
20.7 


Expected Line Shift Associated with Orbital Motion 


| -22 | -322 | -22 | • • • 


Neutral Fc and Ni Fluorescent Lines in GX 301-2 e 


Ni Ka 1.660 
Fc Kf3 1.756 
Fe Ka 1.937 


<f> 0.1-0.2 

nio671 loll 
zlo -610 ls -ll 

25891,0 88% 


4> 0.96-0.98 
221 26 3 5 18 89i 3 43 

k»109 920 24 

285 25 25 788 22 24 


4> 0.47-0.48 

370 17 

-415« 4 60 34« 16 
213 3 _ 8 38 331« 6 





a from House (1969) 

b in units of km s -1 ; positive values arc rcdshiftcd 

c in units of 10~ 5 photons cm _2 s _ 1 

d negative flux value indicates line is seen in absorption 

e for directly observed photons originating on or near the surface of the neutron star 
f GX 301-2 results arc from our study in progress 
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Table 3. Flux and wavelength of H- and He-like lines at 3 phases 



Ion Expected A a 


Ob 1926 


Ob 1928 


Ob 1927 


Ratio 


A 


<t> 


<t> 0.25 




4> 0.5 




Flux 




Shift b Flux c 


Shift 


Flux d 


Shift 


Flux 


1927/26 



H-likc Ions 



sXVI 



Lya 


4.730 


oq 7 342 
oy '-389 


1 2 - 6 






1 90194 

lz,:i -191 


11 1 4 - 5 
il - i -4.6 


9.5 


Si XIV 


















Ly/3 
Lya 


5.217 
6.183 


322 

3341% 


0.3 

1 96.74 
1 - z -0.66 






24275 
z -297 

-163« 42 


4 9 4 1 
18.61% 


16.5 
16.0 


Mg XII 

Ly/3 
Lya 


7.106 
8.422 


333 59 
310«° 59 


1 90.33 
1 - z -0.32 
q fi 0.67 
■ yb -0.51 


1541% 


6 2 3 - 7 
°^-3.5 
5 o,4.6 


44 4 _ 5 64 
_i ko43 
loJ -30 


8 - 8 -2.0 

13-91% 


7.1 
14.4 


No X 


















Ly«5 

Ly7 
Ly/3 
Lya 


9.481 
9.708 
10.239 
12.135 


262 
5 

00 ' -323 

269l 4 47 


0.14 
0.2 

U - 4 -0.36 
1 2 ' 7 


-356 
-186 

cl04 
°-50 


-1.6 

-4.3 2 _1 6 
-1.7 

3-9-2.3 


-132 181 
99C-120 
zz0 -123 

-190l 7 86 
-2641% 


l- 8 -0.7 

4-6l?. 6 

6-3l?.7 
13.7%% 


12.9 
24.2 
15.6 
11.9 


oVIII 


















Lya 


18.970 


563 


0.48 






73 


1.6 


3.4 


He-like Triplets 


sXV 


















r 
i 
f 


5.039 
5.065 
5.102 


401 
426 
264 


0.97 
0.20 
0.44°;t 






iDi -l()8 

-96 

-264«« 4 


8-911.8 

0.6 

4.8 


9.9 
3.0 
10.9 


Si XIII 


















r 
i 
f 


6.648 
6.687 
6.740 


o n7 124 
,3u '-205 

287 
997218 

zz '-293 


1.51°4. 4 
0.1 

0-93% 36 37 


_ 477 275 

4 ' '-293 
152 

-258« 2 65 


-3-9^.9 
6-0ll. 6 

9-0 4 _-3°9 


"184 4 _ 1 124 

-188" 6 94 
-242 5 _ 7 123 


9-41L 
3. 9% 2 
9-4l§. 3 


10.1 

39 

10.6 


Mg XI 
r 
i 
f 


9.169 
9.230 
9.314 


343 120 
°* 3 -152 

309 

354^53 


0.89%™ 
0.22 

0.35°-g. s 


-666^ 

-97 
qoi77 


-5-21L 
4 2 1 - 
4 1 2 ' 3 


_1 7 o46 

-172 

-267 7 _% 


6.01?.! 
6-7l1. 2 


17.3 

8.7 

19.6 


Nc IX 


















r 
i 
f 


13.447 
13.552 
13.699 


201 

299 
328 


0.34 
0.10 
0.42 


359 

153 219 
loo -208 

241 


4.0il 4 


91 ol400 
zlo -830 
-199 

— 941 263 
Z4i -613 


1 1 - 9 

1 6 1 ' 3 
1 - D -1.0 

1 7 - 7 


2.6 

16.1 

5.2 



a from Chandra X-Ray Center Atomic Databasc(http://asc. harvard.edu/atomdb/) 
b in units of km s~~ 1 ; positive values are redshiftcd 
c in units of 10" 5 photons cra _2 s _1 
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d negative flux values indicate lines that are seen in absorption 

Note. — for H-like (and He-like i) doublets the mean A of the 2 lines is used. 
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Table 4. RRC Fit Parameters (redge a model) 



Parameter 


Eclipse 




<t> 0.5 




Nc X 


Nc IX 


Nc X 


Nc IX 


norm b 


0.6 


0.6 


6.5 


6.0 


edge c 


1.3628 


1.1956 


1.3628 


1.1956 


kT d 


0.01035 


0.01035 


0.01035 


0.01035 



a For Nc X the complete fit function used was 
"gauss(l)+gauss(2)+gauss(3)+redge(l)+poly(l)" , and 
for Ne IX "gauss(l)+rcdge(l)+poly(l)"; only redge 
parameters arc displayed 

b in units of 10 -5 photons cm _2 s _1 

c in units of keV 

d in units of keV (kT of 0.01035 keV represents a tem- 
perature of 1.2 X 10 5 X) 

Note. — The measured Lya/RRC flux ratio for Ne X 
was 1.9 in eclipse and 2.1 in <j> 0.5 
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Table 5. Photo-electric Absorption Edges in 0.25 and <j) 0.5 



Edge A(A) Opt. Depth r 



<t> 



0.25 


S K 


5.00 


0.026 




Si K 


6.72 


0.054 




Mg K 


9.50 


0.097 




Ne K 


14.30 


0.42 




Fe L3 a 


17.55 


0.12 


0.5 


S K 


5.0 


0.17 




Fe K b 


1.74 


0.15 



a Fe L3 edge is due to transition from the 
2p 3/2 level 

b Fc K edge is present at 1.74 A in <j> 
0.5, and absent from other phases, consistent 
with variation in the angle of the slab of gas 
with orbital phase. This edge is also present 
(and stronger r = 0.90) in GX 301-2 



